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The effect of pressure on the hexagonal close-packed structure of titanium is investigated. The lattice pa- 
rameters of the equilibrium structure were determined in terms of the Gibbs free energy using the Epitaxial 
Bain Path method. When this process was repeated for several pressures, the effect of pressure on the lattice 

i— i parameters was revealed. The calculated lattice parameters were in good agreement with the experimental 

~—> and theoretical results. The effects of pressure on parameters depending on the electronic structure such as 

*"~5 conductivity and resistivity in the ground state were also investigated up to 30 GPa using density functional 

f— ^ theory. 
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PACS: 61.50.Ks, 05.70.Ce, 72.15.Eb, 71.15Mb 

1. Introduction 



At room temperature and ambient pressure, Ti has a hexagonal close-packed structure called 
the a-phase. The lattice parameters of this structure are a = 2.957 A and c = 4.685 A pQ in 
which the unit cell has two atoms at (1/3, 2/3, 1/4), (2/3, 1/3, 3/4) and the space group number 

Sis 194 (PG^/mmc) with the c/a ratio of ~ 1.59 [2J. Experimental results at room temperature 
indicate that lattice parameters decrease and the c/a ratio increases with pressure [3]. At ambient 

*^J temperature and high pressure, it changes to the w-phase |U|S]. The lattice parameters of this 

structures are a = 4.598 A and c = 2.822 A [1, 6 with three atoms per unit cell at (0, 0, 0), 
(1/3, 2/3, 1/2), (2/3, 1/3, 1/2) and the space group is P6/mmm with the c/a ratio of - 0.61 0. 
The a —¥ u> transition in Titanium is a representative example of martensitic transformation. 
Recently, Trinkel et al. have proposed two pathways for this transformation which are called TAO- 
1 ("Titanium alpha to omega") and TAO-2. This mechanism is a direct mechanism in which six- 
atom transformation proceeds without a meta-stable intermediate phase and has small shuffles and 

^u strains [7]. 

Q^ The effect of pressure on lattice parameters was investigated using the Epitaxial Bain Path 

(EBP) method [HJ E] . In this method, the equilibrium structure is determined by minimized Gibbs 

r- — free energy. As a hexagonal structure is defined by its lattice parameters, the Gibbs free energy 

[G = E(a, c) + PV(a, c)] should be minimized with respect to a and c. The EBP method is 
summarized and explained here. At T = K, the Gibbs free energy was G = E + PV, where P 
is the pressure, V was the volume per atom and E is the internal energy per atom. At P ^ 0, the 
^- lattice parameter of a = a± was chosen and the value of c was varied until one found the values 

c = ci and E = E\ . The energy of E for the lattice parameters c = Cj (i = 1, 2, . . . ) with constant a± 
was calculated using Wien2k until Ci and Ex values satisfied (dE/dc) a = — {Pa\ sin 7) /2, where 7 
is the angle between a and c. When a suitable value of c\ was determined, E\ , V\ and the Gibbs free 
energy G\ , which is G\ = E1+PV1 were found at pressure P. This process was repeated for several 
Oj values at the same pressure and Gi was determined. For the equilibrium structure at pressure 
P, the value giving a minimum Gp was chosen. Thus, by choosing different values of P, the lattice 
parameters a and c, ratio c/a and the volume were determined directly as the functions of pressure. 
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The energy was calculated using Wien2k Package [10], which employs a self-consistent Full- 
Potential Linearized Augmented Plane wave plus local orbital (FLAPW+LO) method, under the 
generalized gradient approximation (GGA) with the Perdew-Burke-Emzerh of 96 exchange Corre- 
lation functional |llj . Moreover, the following parameters were used: Muffin-tin-radius, RMT = 
2.3 bohr, Largest vector in the charge-density Fourier expansion, G max = 12 bohr -1 , K point = 



4954. Plane-wave cutoff. RK„ 



9, cutoff energy = —6 Ryd and charge convergence = 1 x 10 e 



(the charge convergence was used in order to optimize parameters in the SCF cycle). Except for 
K point = 5599 and RMT — 2.247 bohr, other parameters which were related to the effect of 
pressure on the electronic structure did not change. 

2. Results and discussion 

The effect of pressure on the hexagonal close-packed structure of titanium in the ground state 
was investigated. The Gibbs free energy (figure n}, ratio of c/a (figure k|, lattice parameters a and 
c (figure pi and volume (figure B| were calculated as functions of pressure. The experimental data 
were obtained using figure 5 in Errandonea et al. [3J via Gate data software. These results were in 
good agreement with theoretical [5] and experimental studies [3J. 
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Figure 1. The pressure dependence of the Free 
energy G. 



Figure 2. The pressure dependence of ratio 
c/a. 
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Figure 3. The pressure dependence of (a) lattice parameter a and (b) lattice parameter 
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In order to make a comparison between the 
results of this study and those of the refer- 
ences [3], only three points of the pressure (0, 4 
and 7 Gpa) were chosen which was due to the 
limitation of the experimental results in which 
the experimental points were only considered 
up to 14.5 GPa. These comparisons are tabu- 
lated in tables [T] to [3] and shown in figures [5] 
and [6] As shown in the aforementioned figures, 
based on the expectations, the lattice param- 
eters decreased and c/a ratio increased with 
pressure. 

In the present study at T = K and 
P = GPa, a-phase had the c/a ratio of 
1.586 which was in good agreement with the 
experimental 1.584 [3, and theoretical results 
1.611 [S], 1.584 pTJ] or 1.583 [T3] smaller than 
value of 1.633 for an hep crystal. Furthermore, the experimental results showed that the 
for a-Ti gradually increased from 1.584 at atmospheric pressure to 1.622 at 14.5 GPa |3]. 



P(GPa) 

4. The pressure dependence of atomic 



Table 1 . Lattice parameter of hep Ti at zero pressure; estimated deviations with experimental 
work [5] are indicated. 



Lattice parameter 


This work 


Experimental [3] 


Theoretical" [S] 


a (A) 


2.93693±0.007 


2.9575 


2.92455 


c(A) 


4.65834±0.006 


4.68548 


4.71120 


c/a 


1.58613± 0.011 


1.58427 


1.61091 



Table 2. Lattice parameter of hep Ti at 4 GPa; estimated deviations with experimental work [3] 
are indicated. 



Lattice parameter 


This work 


Experimental" [5] 


Theoretical" [8] 


a (A) 


2.89777±0.005 


2.9137 


2.88898 


c(A) 


4.61230±0.008 


4.65088 


4.65735 


c/a 


1.59167±0.003 


1.59621 


1.61211 



Table 3. Lattice parameter of hep Ti at 7 GPa; estimated deviations with experimental work |3j 
are indicated. 



Lattice parameter 


This work 


Experimental [3] 


Theoretical" g] 


a (A) 


2.86973±0.006 


2.88717 


2.86523 


c(A) 


4.58955±0.009 


4.63199 


4.62227 



c/a 



1.59930±0.003 



1.60433 



1.61323 



a Experimental values are obtained by Gate Data Software from figure 5 of the reference [5|. 
b Theoretical values are obtained by Gate Data Software from figure 2 of the reference |Hj. 

The a-phase can transform to the w-phase under pressure. However, it is believed that both 
structures can coexist in the pressure range studied here [S]. The effect of pressure on lattice 
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Figure 5. Lattice parameters vs pressure for 
hep Ti. Estimated deviations with experimen- 
tal work [3] are indicated. 



Figure 6. Ratio of c/a parameters vs pressure 
for hep Ti. Estimated deviations with experi- 
mental work [3] are indicated. 



parameters of a-phase in the range of 0-14.5 GPa was also investigated, both experimentally in |3] 
and theoretically in [8 . These results confirm the possibility of coexistence of both structures in 
the pressure range of 2-9 GPa and even more. 

In fact, the coexistence of these phases has been also reported by experiments within the 
temperature range from room temperature to around 923 K [14_ and by the theory of reconstructive 
phase transitions [15 . 

Table 4. Partial charges in the s, p and d bands. 



P (GPa) 


s 


P 


d 





0.31515 


0.2767 


0.9858 


4 


0.32471 


0.2870 


2.0104 


7 


0.33159 


0.2953 


2.0426 


15 


0.34840 


0.3125 


2.0999 


20 


0.35676 


0.3199 


2.1242 


25 


0.36704 


0.3315 


2.1744 


30 


0.37429 


0.3371 


2.2127 



Table 5. Density of states at the Fermi level 
as a function of pressure. 

P (GPa) n(E F ) (States/ev atom) 




4 
7 
15 
20 
25 
30 



0.8930 
0.8868 
0.8500 
0.8625 
0.8000 
0.8158 
0.7908 



Table 6 

energy. 



. The pressure dependence of the Fermi 



P (GPa) E F (Ryd) 




4 

7 
15 

20 
25 
30 



0.56228 
0.58999 
0.60881 
0.65693 
0.67996 
0.70925 
0.73021 



Using the lattice parameters and Wien2k, the number of electrons in s, p and d bands (tableEl), 
the Fermi energy (table [5| and the density of states at this energy n(e F ) (table [6]) were calculated 
for different pressures. 



23601-4 



Effect of pressure on the electronic structure of hep Titanium 



© 1.25 



£1.15 



P(GPa) 

Figure 7. l/n(e/) as a function of pressure for 
hep Ti. 



Moreover, the electrical conductivity can be 
expressed as a = e 2 TpVpn(eF)/3 |16|[T7] , where 
vf is velocity at the Fermi energy and rp is the 
relaxation time, but ti(£f) has a greater effect 
than the latter two parameters. Figure [7] shows 
1/ti(£f) as a function of pressure. Because p oc 
l/n(ep), figure mean be taken as a measure of 
the effect of pressure on electrical resistivity. 

According to Matthiessen's rule, total elec- 
trical resistivity p, due to electron scattering by 
different factors is given by sum of these factors 
P = Pt\\ + Pd + Pi i where pxh is thermal resis- 
tivity, pd and p\ are resistivity due to defects 
and impurities, respectively. The other factor 
that can scatter electrons is electron-electron 
(e-e) interaction, which is negligible. At high 
temperature, the effects of impurities and de- 
fects are negligible; thus, p sa p Th • a t low temperature, pi-h is less than pn +pi , so p « pu + pi ■ At 
T = K, ions were frozen in fixed positions and electrons do not scatter with phonon (pTh = 0). 
In the present study, pure titanium was investigated. It is chemically and thermodynamically im- 
possible to avoid impurities and defects, so the resistivity at T = is not zero (p = po + Pi 7^ 
where p = p r is called residual resistance). 

Figure |7]shows that at T — K an increase in pressure causes an increase in electrical resistivity, 
which contradicts the results of P.S. Balog, who investigated the phenomenon at 50-700°C, where 
p was due to electron-phonon interaction. A pressure increase (at constant temperature) leads to a 
decrease in inter-atomic spacing and atomic vibrational amplitude, causing a decrease in electrical 
resistivity [13]. However, in the present study, the investigation was carried out at T = K and 
electrical resistivity increased with pressure. According to table HI an increase in pressure leads to 
an increase in the number of electrons per volume. 

Furthermore, according to table [5] an increase in pressure leads to a decrease in the density of 
state, thus causing an increase in the electrical resistivity. 

According to table [6] an increase in pressure leads to an increase in the width of the valence 
band. In the formation of molecules, several atoms are arranged beside each other, so atomic orbitals 
are split and several molecular orbitals are created while the number of orbitals is proportional to 
the number of atoms. However, if the inter-atomic space is small, atomic orbital splitting is larger. 
According to the calculations and figure [2] a pressure increase causes a decrease in the lattice 
constant and inter-atomic space, thus increasing the orbital splitting and width of the valence band. 



Table 7. Partial charges in the p x + p y ,p z ; d z i , d xy + d 7 .i_ i and d xz + d yz 



P (GPa) 


Px + Py 


Pz 


@>xy ~r & x 2_y2 


d z 2 


&xz ~r Uyz 





0.18896 


0.08777 


0.84783 


0.46065 


0.67736 


4 


0.19654 


0.09052 


0.86045 


0.46494 


0.68507 


7 


0.20200 


0.09334 


0.87419 


0.47196 


0.69646 


15 


0.21359 


0.09895 


0.90120 


0.48449 


0.71430 


20 


0.21833 


0.10159 


0.91282 


0.48769 


0.72351 


25 


0.22668 


0.10482 


0.93479 


0.50003 


0.73965 


30 


0.22877 


0.10840 


0.95242 


0.50818 


0.75218 



Table [7] lists the number of electrons in s, p and d orbitals and the deviation from spherical 
symmetry is shown in table [8] This deviation for p and d orbitals is given by fjTJj: 



An d = (n dx 



n d 



{ n d x 



Tic, 



.) 



nd 
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An p = - (n Px +n Py ) - n p __ . 



Table 8. Deviation from spherical symmetry of the p and d states as a function of pressure. 



P (GPa) 


An d 


An p 





0.048 


2.1528 


4 


0.052 


2.1712 


7 


0.054 


2.1641 


15 


0.059 


2.1585 


20 


0.063 


2.1491 


25 


0.064 


2.1625 


30 


0.068 


2.1104 



If An p and And are close to zero, deviation from spherical symmetry will be just slight. Ac- 
cording to table [HJ this deviation increases with pressure for d orbitals. 

3. Conclusion 

The aim of the present study was to investigate the pressure effect on lattice parameters of hep 
structure in titanium. The obtained results showed that the c/a ratio of hep was nearly constant. 
However, it is believed that both structures can coexist in the pressure range studied. The alpha 
phase was the most stable phase at ambient conditions and its transformation to the omega phase 
in the pressure range of 2-9 GPa. Moreover, theoretical and experimental results confirmed the 
possibility of coexistence of both structures within the pressure range of 2-9 GPa and even more. 

Furthermore, effects of pressure on parameters depending on the electronic structure, such as 
conductivity, resistivity, the Fermi energy and n(ep) in the ground state were also investigated up 
to 30 GPa using density functional theory. Moreover, an increase in pressure leads to a decrease in 
the density of state, thus causing an increase in the electrical resistivity. 
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Bn/iMB Ti/iCKy Ha e/ieKTpoHHy CTpyKTypy hep THTainy 

M. fl>Kadpapi, A. /ib+caraH flyer, M. Bae33afle, H. 3apidpi 

TexHO/iorinHi/iw ymBepci/iTeT iM. K.H. Tyci, TerepaH, IpaH 

flocnifl>KeHO BnjiwB Ti/iCKy Ha reKcaroHajibHy mi/ibHO ynaKOBaHy CTpyKTypy TUTaHy. IJapaMeTpw rpaTKH bm- 
3Hanajincn b TepMiHax Bi/ibHoT eHepriT u66ca, BHKopncTOBytoHi/1 MeTOfl eniTaKciT ujjiaxoM BefiHa. Ko.ni/i L4e£n 
npemec noBToptoBaBOR fl/ia fleKi/ibKox twck\b, 6y.no Bi/inB/ieHO Bn.ni/iB Ti/iCKy Ha napawieTpi/i rpaTKi/i. 06hm- 
c/ieHi napaMeTpw rpaTKW flo6pe y3rofl>KyBajiMca i3 TeopeTUHHUMM Ta eKcnepwMeHTa/ibHHMi/i pe3yjibTaTaMM. 
Bn/inB THCKy Ha napaMeTpw, TaKi ak npoBiflHicTb Ta onip b ocHOBHOMy CTaHi, b 3a/ie>KHOCTi Bifl ejieKTpoHHo'f 
CTpyKTypn, TaKO>K 6yno flocjiifl>KeHO a>K flo 30 GPa, BMKopucTOBytoHi/i Teopiro cpyHKL(ioHa/iy rycTMHi/i. 

K/iKJHOBi c/iooa: bujuab Ti/iCKy, hcp-Ti, Bi/ibHa eneprin Ii66ca, napat/ieTpw tpaTKi/t, Merofl EBP 
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